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We examine the reduced phase space of the Bianchi VIIq cosmological model, 
^ ! including the moduli sector. We show that the dynamics of the relevant sector of 
local degrees of freedom is given by a Painleve III equation. We then obtain a zero- 
curvature representation of this Painleve III equation by applying the Belinskii- 
Zakharov method to the Bianchi VIIq model. 
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1. Introduction 

Bianchi VIIq cosmological model is a finite-dimensional dynamical system with 
non-trivial dynamics. It can be considered as a symmetry reduction of the Gowdy 
model [1] whose spatial hypersurface is the three torus T^. The Gowdy models 
represent examples of two commuting Killing reductions of general relativity [2], 
whose space-times have compact spatial hypersurfaces. Although the two commut- 
ing KiUing vector system is integrable [3, 4,5,6,7], it is not known how to extend 
the corresponding methods of constructing solutions [8,9,10,11] to the compact case 
topology. Therefore Bianchi VIIq model can serve as a toy model for exploring the 
issue of integrability of the Gowdy model on the three torus. 

Belinskii and Prancaviglia showed, within a more general framework of Belin- 
skii and Zakharov inverse scattering method [5], that the Einstein equations for 
some Bianchi models admit a zero-curvature representation, indicating that these 
are solvable dynamical systems [12]. However, their analysis was not complete in 
two aspects. The first aspect is related to the fact that their considerations were 
only local, since they ignored topological obstructions coming from the non-trivial 
global topology of the spatial hypersurface, the three torus T^. The problem is that 
spacetimes with compact spatial sections do not allow in general global Bianchi 
metrics, and one can put only locally homogeneous metrics [13,14]. Related to that 
is that compact spatial manifolds have non-trivial topology, and locally diffeomor- 
phic metrics are not necessarily globally diffeomorphic, which means that there are 
global degrees of freedom in the metric, beside the usual local ones. It has been 
shown recently that the moduli parameters enter non-trivially in the diffeomor- 
phism invariant symplectic form, and hence they could change the dynamics of the 
local degrees of freedom [15,16]. The second unexplored aspect is to find out what 

2 



kind of integrable nonlinear dynamical equation can be obtained. 

In order to explore these issues, it is convenient to study the dynamics of 
Bianchi VIIq model in the canonical formalism. We perform constraint and gauge- 
fixing analysis and show that the dynamics of a generic sector of local degrees of 
freedom can be reduced to that of a Painleve III equation. There is also a special 
sector with enhanced symmetry, which has a linear dynamical equation. By using 
the results of Kodama [16] we show that the moduli parameters do not change 
the dynamics of generic local degrees of freedom. We then show how Bianchi VIIq 
model appears in the Belinskii-Zakharov approach, and how to obtain Painleve III 
equation. We then use these results to obtain a linear system whose zero-curvature 
condition is a Painleve III equation. 



Bianchi models are spatially homogeneous spacetimes which admit a three 
dimensional isometry Lie group G that acts simply transitively on each leaf S of 
the homogeneous foliation, for a review and references see [17]. As a consequence, 
there exists for each of these models a set of three left-invariant vector fields L/ on 
E which form the Lie algebra of the group G: 



where C jk are the structure constants of the Lie group. 

Dual to the the vector fields Lj, one can introduce a set of three left-invariant 
one-forms which satisfy the Maurer-Cartan equations 



2. The Class A Bianchi Models 



[Li-, Lj] = C ijLk 



(2.1) 



dx' + i C'jK X^ A x"" = 







(2.2) 
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If the trace C jj of the structure constants is equal to zero, the Bianchi model 
is said to belong to Bianchi class A. For this class of models, the spacetime admits 
foliations by compact slices. 

The structure constants for the class A Bianchi models can always be written 
in the form 

JK = ^jklS^^ , (2.3) 

where ejxL is the totally antisymmetric symbol, and S^^ is a symmetric tensor 
density of weight one over the Lie algebra of G. Further classification of the class 
A Bianchi models is defined with respect to the signature of the symmetric tensor 
density S^"^. The type VII Bianchi model of class A is denoted as Bianchi VIIq, 
and it is characterized by the signature (+,+,0). Hence the structure constants 
for this model are given by 

C^JK = ejifi ^1 + ^JK2 si . (2.4) 

Bianchi models can be considered as the homogeneous sector of general relativ- 
ity. The dynamics can be obtained by performing the corresponding reduction of 
the canonical formulation of general relativity. The canonical variables of general 
relativity are the three metric gij{t,x^) on the spatial section S, and its canoni- 
cally conjugate momenta 7r*-^(i, x*), where are coordinates on E. The action for 
a Bianchi model can be obtained by inserting the expressions 

9ij{t, x) = gij{t)xi{x)xj (x) , n'^it, x) = n'-^{t)L\{x)L^j{x) (2.5) 
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into the canonical form of the Einstein-Hilbert action. This gives 

t 

Sba ^ jdt {i^^Jgij - N^Hi -NHo^ , (2.6) 

to 

where the canonical variables {gjj , tt^"^) have the Poisson brackets 

{ 9IJ : ^''''} - U^f + ^f) ■ (2-7) 

The vector constraint Hj, I = 1, 2, 3, is the reduction of the diffeomorphism 

constraint, and Hj is given by 

Hi = 2 C\i gjL tt^^ - 2 skimS^-^ qjl tt^^ « . (2.8) 

The reduction of the Hamiltonian constraint is Hq, and it is given by 

Hq = (Tr(y7r5r7r) - iTr2(^7r) + T^igSgS) - \Tr:\gS)) ^ . (2.9) 

and N are the corresponding Lagrange multipliers. The constraints (2.8) and 
(2.9) form a closed Poisson algebra 

{Hi,Hj}^C^ijHk , (2.10) 
{Ho,Hi} = , (2.11) 

and therefore they constitute a set of first-class constraints. Furthermore, the equa- 
tion (2.11) implies that the hamiltonian constraint is invariant under the transfor- 
mations generated by the vector constraint. This fact can be used to go to the 
parametrized particle form of the action, defined by the diffeomorphism invariant 
variables and the Hamiltonian constraint. 
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3. Vector Constraint 



The standard approach to finding the dynamics of Bianchi models is to solve 
first the vector constraint. This requires the corresponding gauge- fixing, and be- 
cause of (2.11), one can even find the diffeomorphism invariant variables explicitly. 

In the case of Bianchi VIIq model the topology of E is fixed to be the three 
torus and the coordinates = (x, z) can be chosen such that x^, x^, X" have 
the canonical form 

X""^ = cos 2; dx -|- sin 2; dy , — — sm. z dx -\- cos z dy , — dz . (3.1) 

The three left-invariant one-forms satisfy the Maurer-Cartan equations for the 
type VII Bianchi model 



Note that every Bianchi model has a symmetry group M , which is the auto- 
morphism group of the Lie algebra of G. M can can be realized as a subgroup 
of GL(3 , H) in the following way. Let us consider a set of the left-invariant vec- 
tor fields L/. They form the Lie algebra of G through the commutation relations 
(2.1). An invertible matrix M^j yields a new set of left-invariant vector fields 
Lj = LjMj, whose commutation relations will have the same structure constants 
as Lj if the following identity is satisfied 



dx'+x'Ax^ = 0, dx2+X^Ax' = 0, dx^ = . 



(3.2) 



(3.3) 



where (M ^)^^ is the inverse of the matrix M^j. The matrices Mj with the 
condition (3.3) define the symmetry group M. 
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Equivalently, with the help of the identity (2.3), we can define M as 



S^-' = (detM) ' M\S^^ (M'^') 



1 oKL (T,jT\ J 



L ' 



(3.4) 



where {M'^) l is the transpose of the matrix M'-' i. In particular, in the case of the 
Bianchi-VII model the tensor density S^"^ has the signature (+, + , 0) (see equation 
(2.4)), and hence the condition (3.4) implies the following form of the matrix M 



M = Md Me 




V 








±1 



f COS 9 sin 9 u 
— sin 9 cos 9 v 

V 1 



(3.5) 



Let us consider the following change of variables 



K 



9lj{t)^{ME')j {t)QKL{t){ME)j{t) , 



(3.6) 



where [MeY j{t) is given by 



ME{t) 



( cos9{t) sm9{t) u{t)\ 

-sm9{t) cos9{t) v{t) 
V 1 / 



(3.7) 



and Qij = diag((5i, (^2, Qs) is a diagonal matrix. It is useful to introduce new 
variables (/3o, /?-) as 



<5i = e 



2{/3°+l3++V3l3-) 



2(/3°+/3+-V3/3-) 



2(/3°-2/3+) 



(3.8) 



The definition (3.6) involves only Me, because by rescaling of Q by constants we 
can always put Md — Id. 
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We have two alternative ways of imposing the constraints Hj (2.8). One way 
is to complete the canonical transformation (3.6) by calculating the conjugate 
momenta 

dgij IJ dgij jj dgij jj 

Pe = ^^ ' ^^ = ^^ ' ^^ = ^^ ' ('-^^ 
9giJ IJ dgij J J dgij jj 

and then expressing the constraints (2.8) in terms of the new canonical pairs 
{9 , p0 , u , Pu , V , Pv , (3*^ , Pq , (3~^ , , (3~ , P-). This is the approach taken by 
Kodama [16] . Alternatively, we can calculate tt^"^ as a function oi{9 , u, v , , P^) 
and their time derivatives {9 , ii, v , , /3+ , P~) and then impose the constraints 
Hj. Once the constraints are imposed we can calculate the pre-symplectic form 
and determine the canonically conjugate momenta for the diffeomorphism invari- 
ant variables. We will choose the second alternative, because it is simpler and it 
gives an independent check of Kodama's results. 

Our first step is to express the conjugate momenta n^'^ in terms of the new 
variables 

TT^^ = -Vd^gig'^'KMNg'"' - g" g^'^KMN) , (3.11) 
where the extrinsic curvature Kjj is defined by 

Ku^^{-gij + {Lj^g)ij) . (3.12) 
Our choice of foliation is such that — 0, so that 

Kij^-^gij . (3.13) 



The vector constraint (2.8) is expressed in terms of the product Qirt^^'^ and from 
the equations (3.11) and (3.13) we calculate 

9IK vr^^ = ^ {9IK9^' - 5j {gLK9^')) ■ (3.14) 



Thus, 

Hi^2g2KT^^^^g2Kg^^ , (3.15) 

Hs = 2 (yi;^ 7r^2 _ ^i^i) ^ _ ^2^^^i) . (3.17) 

In order to impose the vector constraint, we only have to calculate the components 
of the product giKg^^ in terms of our new variables (3.6) 

gixg^^ ^ {ME^QMEME-^Q-\ME^y^)^ 

+ [ME^QQ-\ME^y^)\(ME'{ME^Y^)^ . (3.18) 

We substitute (3.7) and (3.8) into (3.18) and after a straightforward calculation 
we find that Hj = is equivalent to 

ii-ve^O , v + ue^O , sinh2(2V^/3~)^ = . (3.19) 
Now we insert (3.6) into (3.11), and by taking into account (3.19), we obtain 
n'^it) = {M^')j'' PklH) mEY)'^j , (3.20) 
where the matrix P = diag(Pi, P2, P3) and it is given by 



P3 = -2^^e-2(/5°-2/3")(/30 + /3+) . (3.21) 

Note that in our calculation of P we have used the constraints (3.19), and therefore 
the expression (3.21) is vahd for all solutions of the constraints. The equations 

(3.19) have two classes of solutions. One class is given by 

u=v=0=Q , 

which corresponds to the sector j3~ 7^ 0, /3~ 7^ 0. This is a generic sector invariant 
under the Bianchi-VII group [16]. In this sector the vector constraint gives that u, 
V and 9 are constants and hence P is automatically diagonal. The second class of 
solutions is 

and it is less obvious that P is diagonal in this case, but it is a consequence of 
the constraints (3.19). Namely, 9 appears in the off-diagonal part of P, but it is 
multiphed by a factor Q^^ — Q2^, which vanishes for /3~ = 0. This is the sector 
which is invariant under a group larger then the Bianchi-VII group [16]. We will 
concentrate on the generic sector, although we will give some brief comments about 
the enhanced symmetry sector. 

The prc-symplectic structure a by definition (2.6) is 

a = 7r^Ugij . (3.22) 

In order to calculate a in terms of the new variables we substitute the expressions 

(3.20) and (3.21) for the three metric gjj and its conjugate momenta n^"^ into the 
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equation (3.22), and we obtain 



a = Tr(Pdg) + 2Tr(PQdM£M£-i) 



(3.23) 



The second term in (3.23) can be calculated from the definition of the matrix Me 
(3.7), and we get 



AMeMe 



-1 



/ d^ (iu-vde\ 

-de d^; + ude 
\ / 



(3.24) 



Prom (3.20) and (3.21) we obtain that the components of the diagonal matrix 



(PQ)'. 



are 



(PQ)\ = 
(PQ)\ = 



\/deig 

N 
Vdet cjf 



(-2/?o + /?+ + V3/?-) , 
(-2/?o + /?■+- a/3 /3i , 



(Pg)\ = -2^^(/30 + /3+) . 



(3.25) 
(3.26) 
(3.27) 



Since PQ is a diagonal matrix, it is obvious that 



Tr(PgdM£M£-^) = . 



(3.28) 



Thus, the second term in (3.23) is identically equal to zero and a is given by 



a = Tr(Pdg) 



(3.29) 



From 



dQi = e2(/?°+/?^+^/^-)2(d/30 + d/3+ + Vsd^) , 



(3.30) 
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dQs = e2 (/^°-2/3+)2(d/30 - 2 d/3+) , (3.32) 



and from (3.21) we obtain 



a ^12 (-/jod/^o + /3+d/3+ + p-dp-^ . (3.33) 

Since 



p„ = - (12 ^)/30 , = (12 ^)p^ . p_ = (12 ^)p- . (3.34) 
the pre-symplectic form becomes 

a = pod/5° + p+d/3+ + p.d/?" . (3.35) 

Note that (3.35) is the pre-symplectic form for the generic sector. The enhanced 
symmetry sector is a special case where p- — (5~ — 0, so that the pre-symplectic 
form becomes 

a = pod/5° + p+d/3+ . (3.36) 

Dynamics of the canonical pairs {(3^ , pQ , (3^ , ps^ , (5~ , p-) is defined by the 
Hamiltonian constraint (2.9). To calculate the expression for the Hamiltonian 
constraint in terms of the diffeomorphism invariant variables we substitute (3.20) 
and (3.21) into (2.9) and obtain 

Ho = (tv (PQPQ) - \ Tr2(PQ) + Tr {SQSQ) - \ Tx^{SQ)^ . (3.37) 

We notice that the Hamiltonian constraint is independent of the unphysical vari- 
ables {9,u,v). In addition, we already have the components of the product 
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{PQYj, which are given by the equations (3.25), (3.26) and (3.27), and SQ is 
a diagonal matrix 

(5g)\ = e2(/3°+/5^+^r), (5g)\ = e2(/3°+/3+-^r)^ {SQ)\ = O . (3.38) 
From this we obtain that the Hamiltonian constraint for the generic sector is given 

by 

^0 = i-m' + + iP-f) + ''^'^ mh^(2x/3/3-) , (3.39) 

or in terms of the canonical pairs (3.34) 



Ho = ^^-^(-(po)' + {p+f + {p-f + 48e4(^°+^^) siTih\2V3(3')) . (3.40) 



The previous calculation applies also to the enhanced symmetry sector, but 
one must take into account that p_ = f3~ = 0. The Hamiltonian constraint then 
becomes 

which corresponds to the dynamics of a two-dimensional relativistic free particle. 
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4. Diffeomorphism invariant phase space 

Note that the reduced pre-symplectic form (3.35) impUes that the difFeomor- 
phism invariant phase space is given by (p^, (3^^) canonical pairs. However, this 
is not correct, since it has been shown recently that the diffeomorphism invariant 
subspace is larger, because it contains the moduli parameters, which are associated 
with the global degrees of freedom of the metric [15, 16]. This is a consequence of 
the fact that Bianchi metrics on compact spatial slices represent compact Reiman- 
nian manifolds with non-trivial topology, and such manifolds can be locally diffeo- 
morphic, but not globally diffeomorphic. 

As a result one cannot define globally a Bianchi metric, so that one needs a 
notion of locally homogeneous spacetime [13,14]. A locally homogeneous Bianchi 
spacetime M with a symmetry group G can be represented as M/K, where M is 
a homogeneous Bianchi space-time with a simply connected spatial section, and K 
is a discrete subgroup of G [15]. In the following we will review Kodama's results 
[16], since his approach is suitable for the Hamiltonian formalism. 

The symmetry group G of the covering space M is defined as 

G = {/ e DiffiM) I - ^} , 

such that it is isomorphic to G, or it is the smallest possible group which contains 
G, and it acts transitively on M. $ = j*$ is the puUback of the canonical data $ 
on M and j : M ^ M is the covering map. 

The moduli space can be parametrized with finitely many parameters {A"}, 
which are determined from the covering maps j\ : M ^ M such that j^(7ri(M)) = 
K\ is isomorphic to K and it is a subgroup of G, where 7ri(M) is the fundamental 
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group of M. One then has to find conjugacy classes of Kx under the homogeneity 
preserving diffeomorphisms HPDG 



HPDG = {/ e Diff{M) I /G/- 1 = G} . 

This procedure gives that K\ = fxKQf^^, where Kq is a reference point, such that 
M/Kq is identified with M, and fx is a finear transformation on M such that 

fix' = FjX' 
flgij = {F^)fgKLF^ 

= {F )kT^ {{F ) )l 

fl^/^g = (detF)(detx)Vdet^7j. (4.1) 
Prom (4.1) one obtains that the pre-symplectic form (3.23) becomes 

a = ^{T^'Jdgu + 2Ca{X)7:"gijd\-) , 
where V — Jj^^ (fix det x, 

n{X) = J d^x{detx){detF) , Ca{X) = ^ J d^x{detx){det F)daFF-^ , 

Do Do 

and Dq is the fundamental region of the action of Kq on M. The Hamiltonian 
constraint is rescaled by V/il{X). 

In the case of Bianchi VII model, G — VIIq, where + fixes the orientation of 
the spatial section. M is a three-torus T^, which can be represented as E^/K where 
is the Euclidian space, and K — Z^. G — GxD2, where D2 is the dihedral group 
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and X denotes a semi-direct product. K\ are represented by GL{'i^'R) matrices, 
whose last row is given by (Ztt, rnvr, nvr), l,m,n & Z. By using HPDG and modular 
transformations, can be put into form 



/X Y 2 

W \ , X>0 , n>0 

\ nn, 



(4.2) 



Therefore X, F, Z, W and n represent the modular parameters. The Kq is associ- 
ated with X — 1,Y — Z — W — 0, and the deformation map fx is given by the 
matrix 



/X Y Z\ 

x-^ w 
Vo 1 / 

The fundamental region is Dq — {0 < x,y < 1,0 < z < nn}. 
From (4.3) one obtains 



(4.3) 



F = R3iz)AR3{-z) , 



(4.4) 



where Rs is the rotational matrix around z-axes. The correction in the pre- 
symplectic form due to modular parameters is given by 

Tr{FF-^pg) = cos(2^) - \{XY - YX) sin(2^)] (QiPi - Q2P2) ■ (4.5) 

However, when (4.5) is integrated over Dq, one obtains zero, and hence 5a vanishes 
[16]. Therefore the moduli parameters do not enter into the reduced pre-symplectic 
form, and hence they will not influence the dynamics of the local degrees of freedom. 
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5. Hamiltonian Constraint 

The action for the generic sector of the Bianchi VIIq model can now be ex- 
pressed in terms of the canonical variables {pnil3^) as 

t 



were N is the Lagrange multiplier and — V^et qEq/A is the rescaled Hamil- 
tonian constraint (3.40). This system is reparametrization invariant i.e. it has a 
symmetry generated by the constraint Hq 



where e is the parameter of the transformation. 

In order to find the dynamics of the physical degrees of freedom, we need to 
fix the reparametrization gauge symmetry (5.2). As discussed in [18], this type of 
gauge-fixing requires the specification of the time variable, in addition to the usual 
requirement that the Faddccv-Popov determinant is non-zero. In the case of our 
system, the analysis simplifies if we introduce new canonical coordinates 




(5.1) 



dN 



(5.2) 



T 



4 (/3" + /3+) 



6(/?o - /3+) , 



In 16, PT = - (po+P+) , 

PI = ^ (po - P+) , 
1 



(5.3) 




The Hamiltonian constraint now becomes 



■PTPi + + 8e^ sinh2(-) , 



(5.4) 
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Since (5.4) is independent of q , it follows that pi is a constant c. On the other 
hand, the Lagrange multiplier N plays the role of the one-dimensional metric on 
the world-line parametrized by t, and hence it can be set to a positive constant via 
(5.2). Since the equation for T is given by 

f + cN^O , (5.5) 

then one can choose the following gauge 

T^t , N^-l/c . (5.6) 

Note that the gauge choice (5.6) requires pi = c < 0, in order for to be positive 
and finite. Configurations with pi — —c are physically equivalent to pi = c config- 
urations, because our system is invariant under time-reversal (configurations with 
pi = —c have T = —t). The configurations with pi = are excluded, since they 
belong to the enhanced-symmetry sector p- — (5~ — 0. 

The equation (5.6) defines the required gauge choice, and therefore the Hamil- 
tonian of the physical degrees of freedom is given by solving the Hamiltonian 
constraint for pj- 



PT 




Hence the physical phase space is given by the (p, q) canonical pairs, and the 
corresponding Hamiltonian is given by H* . The dynamical consistency of the 
gauge choice (5.6) can be checked explicitly, by comparing the equations of motion 
for pt-iP-i q from (5.4) with the corresponding equations coming from the reduced 
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Hamiltonian H*. For example 

PT = -8A^e smh (-) = + {H ,pt} = , (5.8) 

where {, }* is the Poisson bracket with respect to the reduced phase space variables 

Prom H* we get p — q and p — — 4e* sinh q, so that 

(fa 

+ 4e* sinh q = . (5.9) 

This is a Painleve III equation [22]. One can put (5.9) into the standard form via 
time redefinition t — 

d / dq. ^ . , 
— (r— 1 + 4smhg = 
or ^ or ^ 

Note that the standard form (5.10) could have been also obtained directly by 
choosing the gauge T = logt, pi = —c and N = —l/{ct). 

6. Zero-curvature representation of Painleve III 

By using a more general framework of Belinskii-Zakharov inverse scattering 
method for the spacetimes admitting two commuting spacelike Killing vectors [5] , 
Belinskii and Prancaviglia showed that the Einstein equations for certain Bianchi 
models admit a zero-curvature representation [12]. In this section we will review 
their approach for Bianchi I, II, VIq and VIIq models. A particular attention will 
be given to the Bianchi VIIq model. We will show that within the framework of the 
inverse scattering method, the dynamics of the Bianchi VIIq model is given by the 
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(5.10) 



Painleve III equation (5.10). This is in complete accordance with the results from 
the previous section. In addition, we will derive a zero-curvature representation 
for this Painleve III equation. 

We begin with a brief discussion of Belinskii and Zakharov method for the 
midi-superspace models that are characterized by the existence of a two-parameter 
Abelian group of motions with two spacelike Killing vectors [5] . Let us choose co- 
ordinates adapted to the action of the symmetry group so that the metric assumes 
the following form 

^ _f ^^2 ^^2^ ^^a^^b ^ ^Q -^^ 

where a,b = 1,2, {x^ , , x'^ , x'^} = {t,x,y,z}, / is a positive function and gab 
is a symmetric two-by-two matrix. The function / and the matrix g^i, depend 
only on the co-ordinates {t,z}, or equivalently on the null co-ordinates {C,v} — 
{^{z + t),^{z-t)}. There is a freedom to perform the co-ordinate transformations 

icM^mim}- (6.2) 

It is easy to see that the transformations (6.2) preserve both the conformally flat 
two-metric f{—dt^ -\- dz^) and the positivity of the function / if d^^ drjfj > 0. 

The complete set of vacuum Einstein equations for the metric (6.1) decomposes 
into two groups of equations [5]. The first group determines the matrix g^b and 
can be written as a single matrix equation 

drj {a d^g g~^) + {a dr^g g~^) = , (6.3) 
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where c? = det g and rj} are the null co-ordinates. The second group of equa- 
tions determines the function f{^,r]) in terms of a given solution of (6.3): 

df(\na) 1 _ 
9^(ln/) = + -^trA\ (6.4) 

9^(lna) 1 o , . 

9, In/ = /—^ + - tvB^ , 6.5 

where a^ — d^a, arj — dr^a and the matrices A and B are defined by 

A^-a d^g g~^, B ^ a dy^g g~^ . (6.6) 

The dynamics of the system is thus essentially determined by the equation (6.3). 
By taking the trace of the equation (6.3) and by using the definition for a, we 
obtain 

a^fj = . (6-7) 
The two independent solutions of this equation are 

a = c(0 + d{r)) , P = c(0 - d{r]) . (6.8) 

By using the transformations (6.2), one can bring the functions c(^) and d{r]) to 
a prescribed form. However, we will consider the general form without specifying 
the functions c{^) and d{r]) in advance. 

The crucial step in the inverse scattering method is to define the linearized 
system whose integr ability conditions are the equations of interest, in our case the 
equation (6.3). Following ref. [5], we define the two differential operators 

Di = d^ -^dx , (6.9) 
A — a 
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D2^d, + ^ dx , (6.10) 

where A is a complex parameter independent of the co-ordinates {C,v}- K is 
straightforward to see that the differential operators Di and D2 commute since a 
satisfies the wave equation (6.7) 

[Di, D2] = a^^ dx^O. (6.11) 

The next step is to consider the following hnear system 

Dii/j^-^jjj, (6.12) 



A — a 

D2^ = T^^, (6-13) 



where ilj{\,^,ri) is a complex matrix function, and the real matrices A, B and 
the real function a do not depend on the complex parameter A. The integrabil- 
ity conditions for the system (6.12) and (6.13) are given by the equation (6.3). 
Furthermore, a solution 'ip{X,^,7]) yields a matrix g{^,ri) that satisfies the original 
equation (6.3). Namely, the matrix g{^,r]) is given by 

9{tv)-^{\^:V)\x^o- (6.14) 

In order to take into account that g{C,v) is real and symmetric we have to impose 
two additional conditions, see [5]. Also, it is easy to see that the equations (6.12) 
and (6.13) for A = 0, imply equations (6.6). 

Although Belinskii and Prancaviglia formulation is more general [12], we will 
discuss only type A Bianchi models which are compatible with the inverse scatter- 
ing method. It is not difficult to show that Bianchi types I, II, VIq and VIIq admit 
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the representation (6.1). Namely, the metric for these Bianchi types has the form 



ds^ = -dT^ + g.. dx^dx^ , (6.15) 



where 



9ij = 9lJ X iX j ■ (6.16) 

For these models it is always possible to have the one forms in the following 
form 

= l\dx + l\dy , x^ ^l^idx + f2dy , X^ = dz , (6.17) 
where I'^i, are functions of z only. Let us consider the two-by-two matrix 



1 f 2 



(6.18) 



An important consequence of the Maurer-Cartan equations for the one forms x^ is 
that the matrix I satisfies the following linear differential equation 

^^C^d , (6.19) 
dz 

where the matrix C is the same matrix as the upper two-by-two block on the princi- 
pal diagonal of the matrix S^"^ defined in the equation (2.3). e is the antisymmetric 
matrix with ei2 — 1- 
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After a time redefinition t = t{T), the metric (6.15) can be written in the form 
ds"^ = f{t) {-dt^ + dz^) + gab{t, z) dx^'dx^ . (6.20) 
Here / is a function of t only, and 

g{t,z) = f{z)^{t)l{z) , (6.21) 

where I is given by (6.18) and 7 is a two-by-two symmetric matrix. Notice that 
now 

= (det0^det7 . (6.22) 

Moreover, for these models, the determinant of the matrix I is always equal to one, 
i.e. det Z = 1, so that 

a'^{t) = det7(t) . (6.23) 
In addition, a has to satisfy the equation (6.7), which now reads 

a(t) = . (6.24) 

Hence, a can only be a linear function of time. 

As Behnskii and Prancavigha have showed [12], the hnearized system (6.12) 
and (6.13) can be simpHfied for the models described by the metric (6.20). The 
first step is to define a two-by-two matrix function ip by 

jjj^fifil, (6.25) 
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and a constant two-by-two matrix 



R^eC . (6.26) 

The second step is to substitute (6.21) into (6.6) and use the definition of the coor- 
dinates ^ and 77. Then the results of these calculations, together with the definition 
(6.25), can be used to simplify the equations (6.12) and (6.13). The crucial step in 
which a simplification occurs is the coordinate transformation {i, A} — > {t, A}, 
where w is given by 

2 

w = l{\ + 2l3 + \) . (6.27) 
A 

To perform this co-ordinate transformation we can use a and (5 as given by (6.8). 
The linear system after this co-ordinate transformation involves only derivatives 
in t and A since all the terms involving derivatives in w are canceled. Finally, it 
is useful to make some simple linear combinations of the two equations and to use 
the fact that a is a linear function of time. In this way one obtains a new linear 
system 

dx^ = ^ {-R^ - + ^77" V + -^^R^ - ^7i?^7" V) • (6.28) 

Although the matrix function (p{t,\w) depends on all three variables, the right- 
hand side of the system (6.28) does not have any w dependence. 

The integrability condition for the system (6.28) is 

--^(a77~^) = R-fR^-f-^ - 'fR^-f-^R . (6.29) 
a at 

25 



To derive the equation (6.29) from the system (6.28) it is necessary to use the fact 
that for these models a is a hnear function of time. 

Equivalently, one can derive the equation (6.29) by a direct substitution of 
the formula (6.21) into equation (6.3). A straightforward calculation, using the 
definition of ^, rj, the equation (6.19) and the fact that a is a function of time 
only, yields the equation (6.29). Thus the dynamics of the these Bianchi models is 
essentially determined by the equation (6.29). Furthermore we have confirmed that 
the linear system (6.28) corresponds to the Bianchi models under consideration. 

Let us now consider Bianchi VIIq model. In that case the spatial hyper-surface 
is a three torus T^. As we have shown, the modular parameters do not affect the 
dynamics of the local degrees of freedom, and hence the equations (6.29) and (6.28) 
will be correct dynamical equations. 

The matrix I for Bianchi VIIq model is given by 

(cos z sin z 
— sm z cos z 

and therefore R is 




(6.31) 



We can take 7 to be diagonal 




(6.32) 



We also choose a = t, thus the following relation between the functions a and h 

a — ah — t . (6.33) 
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In order to derive the differential equation which defines the dynamics for this 
model we substitute formulas (6.31) and (6.32) into equation (6.29) and use the 
relation (6.33) in order to eliminate the function b. A straightforward calculation 
yields the following scalar equation 

-^{2td/a) = t'^a-^ - aH-^ . (6.34) 

The redefinitions r = t^/A and = a'^/Ar then give the Painleve III equation 
(5.10). If we define / = then the equation (5.10) becomes 



The equation (6.35) is the canonical form of the Painleve III equation with the 

coefficients a = —2, (3 = 2, ^ = 6 = 0, see [22]. 



The dynamics of the generic sector of the Bianchi VIIq model is given by a 
Painleve III equation, and therefore it is an integrable model. The moduli pa- 
rameters do not affect the dynamics of this sector, and hence do not spoil the 
integrability. Therefore the original claim by Belinskii and Francaviglia that the 
Bianchi VIIq model is integrable is shown to be correct. Furthermore, by using 
their results, we have found a zero- curvature representation of the corresponding 
Painleve III equation. 

In the enhanced symmetry sector we have obtained a linear dynamical equation. 
We expect that the moduli parameters would not spoil this, although a thorough 
investigation of this point would be necessary. 




(6.35) 



7. Conclusions 
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The result that a Painleve III equation appears as the dynamical equation 
of the local degrees of freedom in the Bianchi VIIq model can be used to obtain 
information both about the Painleve III equation and about the physical properties 
of the model. 

As far as the theory of Painleve III equation is considered, the linear system 
(6.28) represents a new tool for the study of Painleve III equation. This linear 
system is different from the linear system which is used for the study of Painleve 
III equation within the isomonodromic deformation method [19,20,21]. However, 
it remains to be explored what are the advantages of the new linear system. 

On the cosmology side, one can now examine the physical properties of the 
solutions, like small and large time asymptotic, as well as the singularities, since 
these properties of the Painleve III solutions have been thoroughly studied [21]. 
In addition, the quantization of this model should be straightforward, since the 
reduced phase space Hamiltonian can be promoted into a Hermitian operator. 
However, it is not clear whether one can find exact solutions for the quantum 
dynamics, since the Hamiltonians at different times do not commute. 
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